A hybrid yang transform adomian decomposition method for solving time-fractional nonlinear partial differential equation

Nonlinear time-fractional partial differential equations (NTFPDEs) play a great role in the mathematical modeling of real-world phenomena like traffic models, the design of earthquakes, fractional stochastic systems, diffusion processes, and control processing. Solving such problems is reasonably challenging, and the nonlinear part and fractional operator make them more problematic. Thus, developing suitable numerical methods is an active area of research. In this paper, we develop a new numerical method called Yang transform Adomian decomposition method (YTADM) by mixing the Yang transform and the Adomian decomposition method for solving NTFPDEs. The derivative of the problem is considered in sense of Caputo fractional order. The stability and convergence of the developed method are discussed in the Banach space sense. The effectiveness, validity, and practicability of the method are demonstrated by solving four examples of NTFPEs. The findings suggest that the proposed method gives a better solution than other compared numerical methods. Additionally, the proposed scheme achieves an accurate solution with a few numbers of iteration, and thus the method is suitable for handling a wide class of NTFPDEs arising in the application of nonlinear phenomena.

In most cases, fractional differential equations (FDEs) are considered a generalization of differential equations (DEs) because they describe functional values at every continuous point without losing any memory or hereditary behaviors of natural phenomena.For this reason, many authors have studied FDEs for modeling and deep understanding of real-world natural phenomena such as business models with fractional derivations, evaluation of viscoelastic surface characteristics, human behavior in the field of mathematical psychology [27,28], models of love between humans [29], and models of happiness and in other fields of applications [23-25, 30, 31].
Nonlinear time-fractional partial differential equations (NTFPDEs) have played an alternative role in diverse fields such as applied science, mathematics, physics, medicine, biology, geology, chemistry, and engineering due to the non-locality nature of the power-law structure with arbitrary positive order.In the modern age, it is impossible to model many real-world problems without using nonlinear FPDEs.Fractional calculus can be called this century's calculus because of the diversity of applications in every discipline of applied science, technology, and engineering [32][33][34].
In applied mathematics and many other areas of science, it is important to obtain accurate solutions to FPDEs.In general, FPDEs do not have closed-form exact solutions for most problems, and it is very difficult to obtain exact solutions [35].For this reason, various numerical methods have been developed for solving linear and nonlinear FPDEs.For instance, Mohamed and Torky [33] proposed the Legendre multiwavelet Galerkin method to solve the conformable fractional modified KdV-ZK equation.Singh et al. [36] announced the homotopy perturbation Sumudu transform method for nonlinear FPDEs arising in the spatial diffusion of biological populations.Wang and Liu [37] familiarized themselves with the new Sumudu transform iterative method for the time-fractional Cauchy reaction-diffusion equation.Ziane et al. [38] introduced the Elzaki transform combined with the variational iteration method for PDEs of fractional order.Dehestani et al. [39] proposed Legendre-Laguerre functions based on Legendre and Laguerre polynomials for solving a class of linear and nonlinear time-space fractional PDEs with variable coefficients.Dehestani et al. [40] presents the discrete scheme based on Genocchi polynomials and fractional Laguerre functions to solve multiterm variable-order time-fractional PDEs in the large interval.Dehestani et al. [41] familiarized themselves with the fractionalorder Genocchi-Petrov-Galerkin method to investigate the approximate solution of time-space fractional Fokker-Planck equations.Wang et al. [42] proposed two analytical methods, the residual power series method and the homotopy analysis transform method, to solve NTF-PDEs with proportional delay.
Recently, researchers have used several numerical methods to solve different types of NTFPDEs.For example, Bekela et al. [32] applied the hybrid numerical method of Laplace-like transform and variational theory for solving NTFPDEs with proportional delay.Malyk et al. [43] employed the Yang-Abdel-Cattani derivative operator for solving the nonlinear fractional diffusion equation.The numerical method is presented by using a novel spline technique for solving fourth-order timefractional evolution problems [44].Al-Deiakeh et al. [45] introduced an approximation for nonlinear FPDE using the combination of the multi-Laplace transform and ADM.Mallick et al. [46] proposed an iterative method for solving time-fractional PDEs with proportional delays.A wide application of NTFPDEs in many real life phenomena and the growing of to search for suitable numerical methods motivate us to propose a hybrid numerical method called the Yang transform Adomian decomposition method (YTADM) for solving NTFPDEs.
George Adomian introduced a modern mathematical method to solve nonlinear DEs in the 1980s, described as ADM [47][48][49].Similarly, another powerful method found by Yang to solve linear PDEs was described as YT, which transforms the initial DEs into an algebraic equation [50][51][52].However, it does not handle nonlinear terms in NTFPDEs because of the difficulties that are caused by nonlinear terms.ADM is a powerful numerical method for solving nonlinear DEs, but it does not properly decompose the fractional order parts.The main objective of this study is to combine two methods for solving the NTFPDEs: YT is used to decompose FDEs into algebraic equations, and ADM is used to decompose the nonlinear terms in NTFPDEs into a series solution.The stability and convergence conditions of the proposed method were analyzed.In addition, some test examples of NTF-PDEs are solved to illustrate the capability, accuracy, and simplicity of the proposed method.
Furthermore, to demonstrate the versatility and robustness of the YTADM in different areas of science, we take four examples that have a wide range of applications in various fields of science and technology.Accordingly, first we considered the nonlinear time-fractional Newell-Whitehead-Segel equation (NWSE), which predicts the existence of traveling wave patterns.These kinds of systems can be observed in a variety of natural systems, including physical, biological, and chemical systems.For example, in medicine: the mechanisms of spreading epidemic diseases, healing wounded tissue, chemical reactions of multiple substances, and neutron diffusion theory all belong to reaction-diffusion systems [49].
Then, we considered the nonlinear fractional Cauchy reaction-diffusion equation (CRDE), which is used to describe the evolution of a system over time with reaction and diffusion processes in various fields.For instance, in biology and medicine, these equations are used to model the spreading of diseases, tumor growth, etc.In chemistry, CRDEs are employed to understand chemical reactions that exhibit complex kinetics and diffusion processes.In population dynamics, these equations are employed to model the behavior of ecological and natural populations.In material science, CRDEs are expended for modeling transport marvels in porous media whereas in finance and economics, they are used for modeling the financial markets with memory goods and long-range dependence [47,48], and references therein.
Hereafter, we solve the nonlinear time-fractional Fornberg-Whitham equation (FWE).The FWE is a type of traveling wave equation and have variety of applications in physics and engineering.Such as in the propagation of electrical signals and optimization of guided communication systems.In mathematical physics to investigate how non-linear dispersive water waves break [50].Final, we considered one problem from nonlinear time-fractional Fokker-Planck equation (FPE).The nonlinear time-fractional FPE arises in various fields like in chemistry, natural science, including astrophysics problems, biological applications, chemical physics, and other fields [51].
Here, the numerical results obtained by YTADM are summarized in tables and figures to illustrate easily understanding for various fractional order α .For the fractional order α = 1 , the obtained numerical solution have excellent agreement with the exact solution.It is evident from the theory of fractional calculus that the approximate solution of the problem with derivative 1 continuously tends to the exact solution when the fractional derivative α(m − 1 < α ≤ m) tends to positive integer m = 1 [52].
The remaining paper is organized in the following way: Sect."Preliminaries" of this work gives detailed concepts related to fractional integrals, derivatives, and YT.Sect."A hybrid numerical method for solving NTFPDEs" familiarizes reader with the mathematical formulation of YTADM using YT and ADM.Sects."Stability analysis of YTADM to solve NTFPDEs" and "Convergence analysis of YTADM for solving NTFPDEs" offer stability and convergence of the proposed method in the Banach space sense.Sect."Numerical results and discussions" presents a numerical simulation of the method on some NTFP-DEs.Finally, the conclusion of the method is presented in Sect."Conclusion".

Preliminaries
This part is devoted to some basic concepts and definitions of YT, fractional integrals, and fractional derivatives, which are essential for accepting the remainder of the monograph.Definition 2.1 [6,7].A real function f (t), t > 0 is in the space C τ , τ ∈ R if there exists a real number (p > τ) , such that f (t) = t p f 1 (t) , where f 1 (t) ∈ C[0, ∞) , and it is said to be in the space Definition 2.2 [24,25].The left sides Riemann-Liouville fractional integral of order α ≥ 0 , of a function where Ŵ(.) is the gamma function.Definition 2.3 [7,55].The fractional derivative of f (t) in the Caputo sense is defined as.
where a C D α t is the Caputo fractional derivative operator of order α , n − 1 < α ≤ n, n ∈ N. Definition 2.4 [56,57].The YT is defined over the set of function.
where C is constant and m 1 , m 2 are either finite or infinite by the following integral formula It converges if the limit of the integral exists, and diverges if not.Accordingly, the following results follows directly from the above integral formula of YT [56,59]: linearity property.Definition 2.5 [56,59].The YT of the n th derivative of the function f (t) is defined as Definition 2.6 [58].The YT of the function f (t) in the Caputo fractional derivative is defined as. for Definition 2.7 ([32]) Let P : X → X be a mapping of a set X into itself.The fixed point of P is x ∈ X which is mapped onto itself, that is, the image Px coincides with x.Definition 2.8 (Contraction [32]) Let X = (X, d) be a metric space.A mapping P : X → X is called a contrac- tion on X if there is 0 ≤ ε < 1 such that, for all x, y ∈ X, Theorem 2.1 (Banach's fixed point theorem [32]) Consider a metric space X = (X, d), where x = ∅.Suppose that X is complete, and let P : X → X be a contraction on X.Then P has a unique fixed point.Theorem 2.2 [60].Let (X, .) be a Banach space and P be a self-map of X (P : X → X), satisfying.
then P is Picard, P-stable.

A hybrid numerical method for solving NTFPDEs
In this section, we derive a hybrid numerical method that uses YT and ADM.We named this developed method YTADM.To illustrate the basic idea of this method, we consider the following general NTFPDEs with a source term [38,58] with the initial condition where 0 C D α t = ∂ α ∂t α is the Caputo fractional derivative of order α with respect to t , R is the linear differential opera- tor, N represents the general nonlinear differential opera- tor and f (x, t) is the source term.To develop YTADM for solving (1) and ( 2), we follow the following steps.
Step 1: Apply YT on both sides of (1) to obtain By using Definition 2.6 we transform the fractional derivative on the left side of (3) in to Note that in this step, YT changes the fractional derivative into an algebraic equation, and using the given initial condition, we can rewrite (4) as where By the linearity property of the YT, the above equation becomes Step 2: Apply the inverse YT on both sides of (5), to obtain where Step 3: In order to decompose the nonlinear term in NTFPDE and generate a series solution of the proposed method, the ADM is applied.The ADM defines the solution in an infinite series of linear terms as and the nonlinear term in the problem expressed in terms of the Adomian polynomial as follows: where Step 4: Substituting the obtained values in ( 7) and ( 8) into ( 6), we obtain Step 5: To obtain the component solutions of the proposed method, we equating the terms on both sides of (10), and get the following relation Step 6: Following component approximate solution, we obtain the general solution in a series form as Note that it is impossible to find infinite component of solution.Hence, after (m + 1) th truncation, the iteration formula of YTADM is given by ( 6) which is the (m + 1) th approximate solution of the NTF- PDE in (1) for 0 < α ≤ 1.

Stability analysis of YTADM to solve NTFPDEs
In this section, we discuss the stability analysis of the method presented in Sect."A hybrid numerical method for solving NTFPDEs".For this, we state and prove an essential condition concerning the stability of YTADM when it is applied to solve NTFPDEs.To reveal Picard stability, it enough to show that the mapping associated with YTADM in (13) satisfies the conditions of Theorem 2.2.
Theorem 3.1 Let (X, .) be a Banach space and P : X → X be a self-map of X .Then, the iteration proce- dure of YTADM defined by.
is P-stable if.
Proof.First we show that P has a fixed point.To do this, for n, m ∈ N , we have.
By subtracting (15) from ( 14), we obtain ( 14) Taking the norm on both sides of ( 16), without loss of generality, we have Using the linearity property of the YT and its inverse, we obtain.
Using the properties of the norm, we proceed.Now, assuming that.
Then, (17) becomes However, from the properties of YT, we have.
Therefore, from (18), we obtain where ε = (ε 0 + ε 1 ) . Hence, the self-mapping P has a fixed point.Now, we show P satisfies the condition in Theorem 3.1 .For this, we have This shows that the conditions of Theorem 3.1 hold for self-mapping P .Hence, by Theorem ,

Convergence analysis of YTADM for solving NTFPDEs
In this section, we prove the convergence of YTADM when applied to NTFPDEs.
Theorem 3.2.Let (X, .) be a Banach space and P : X → X be a mapping associated with YTADM be defined by (13).Then, P has a unique fixed point and the sequence solution Proof.First, we show the existence of a fixed point of P .For this define that {s n } is the sequence of partial sums of the series Eq. ( 13) as, Now, we show that {s n } ∞ n=0 is a Cauchy sequence in Banach space (X, .) .For this purpose, we consider.
Now, for every m, n ∈ N with n > m , then by using (20) and triangle inequality, we obtain Since 0 < ε < 1 , the sum 1 n=0 is a Cauchy sequence and is convergent.Hence, P has a fixed point.
Next, to complete this proof, let {u n (x, t)} converge to u ∈ X and we need to show that u is the unique fixed point of P .To achieve this, let l be another fixed point of P .Then, by (19), we have Since,(1 − ε) < 0 for 0 < ε < 1 , the above inequality can only hold if �u − l� = 0 implies u = l .Hence, u is the unique fixed point of P , which completes the proof.

Numerical results and discussions
In this section, four examples of NTFPDEs are solved by using the YTADM.The validity of the proposed numerical method is observed in terms of its absolute errors.To show the behaviors of the corresponding examples, some figures are also plotted for fractional orders α .All the results are calculated using MATLAB.

with initial conditions
The exact solution of this problem is for the special case α = 1 .For sim- plicity, we drop (x, t) from u(x, t) and ∂ 2 ∂x 2 u(x, t).
Step 1: Applying YT on both sides of (23) and using the differentiation property, we have In our case 0 < α ≤ 1 , ( 24) is reduced to Using the given initial condition, we obtain Step 2: Applying the inverse YT on both sides of (25), we obtain.This implies that.
Step 3: ADM defines the solution by an infinite series of linear terms as and the nonlinear term as where A n (u) are the Adomian polynomials to be determined.
Step 4: substituting ( 28) and ( 27) into (26) by using the linearity property of the YT, to obtain u(x, 0) = . (24) Step 5: Equating the terms on both sides of (29), we have the following recurrence relation.By using ( 9), the few components of the Adomian polynomials for the nonlinear term u 2 have been derived in the form Consequently, solving the above equations, the first few components of the YTADM solution for (30) are derived as follows: Using the same technique, , and so on.
As shown in Figs.1a and b, the approximate solutions of the pattern behavior gradually decrease as the fractional order values get closer and closer to.The same as in the first example, for comparison purposes, we consider only the special case because the exact solution exists.For this reason, in Fig. 2, we compare the approximate solution with the exact solution.
As seen from Fig. 2, the approximate solution obtained is similar to the exact solution.This finding implies that the results obtained by YTADM are in good agreement with the exact results.The numerical solution to this problem is also compared with the results of the Laplace transform decomposition method (LTDM) [64], the Aboodh transform homotopy perturbation method (ATHPM) [65], and the q-homotopy Shehu analysis transform method (q-HSATM) [61] in Table 1 at the same number of iterations and domain discretization.
Table 1 shows the numerical approximate and comparison solutions of YTADM and other numerical methods.The absolute errors of YTADM at various points in the corresponding domain are presented and compared with LTDM [64], q-HSATM [61], ATHPM [65], and the exact solution at α = 1 , it yields far more (31) u(x, t) = u 0 (x, t) + u 1 (x, t) + u 2 (x, t) + . . ., robust results than all the other numerical methods.Hence, we conclude that the absolute error determined in the table shows that the method is too accurate for treating NTFPDEs.

with initial condition
The exact solution of this problem is u(x, t) = e x+t for the special case α = 1 .For simplicity we drop (x, t) from each term of (32).To solve this problem by YTADM, we follow the procedure stated in Sect."A hybrid numerical method for solving NTFPDEs".
Step 1: Apply YT on both sides of ( 32) to obtain In our case 0 < α ≤ 1 , ( 33) is reduced to Step 2: Taking the inverse YT on both sides of ( 34), ( 32) Step 3: The ADM defines the solution by an infinite series of linear terms as and the nonlinear term as (36) where A n and B n are the Adomian polynomials to be determined.
Step 4: Substituting ( 36), ( 37) and ( 38) into ( 35) by using the linearity property of the YT, to obtain (38)   Step 5: To obtain the component solutions, we equating the terms on both sides of (39) to get By using (9), the few components of the Adomian polynomials for the nonlinear terms uu xx and u 2 have been derived in the form and so on.
Step 6: Following component approximate solution, we obtain the general solution as where E α (.) is well-known as the Mittag-Leffler function defined in [7].
To show the solution behavior of the approximate solution obtained by YTADM, we consider different values of the fractional order α (α = 0.4, 0.6, 0.8, 1) and depict the results in Figs. 3 and 4.
The numerical solutions obtained by YTADM for distinct fractional orders are depicted in Figs.3a and (39) We have observed from this figure that the solution behaviors of the time-fractional derivatives are clearly decreasing as the values of fractional order increase.In the case of ordinary derivatives, the solution behavior of such a DE is not clearly shown.Therefore, arbitrariness in fractional-order derivatives introduces more degrees of freedom in the design and study of realtime events.It is observed from Figs. 4b and d that the approximate solution obtained by YTADM is almost identical to the exact solution at.From Fig. 4, it is found that exact and approximate solutions are in complete agreement.The absolute errors illustrated in Figs.4c and e indicate that the designed numerical method is a suitable one for solving nonlinear fractional Cauchy reaction-diffusion equations that arise in various fields of science.We have also compared YTADM with the Aboodh variational iteration method (AVIM) [63] using their absolute errors.For this reason, we consider the time discretization points as in AVIM and record the obtained results in Table 2. Table 2 shows the numerical solution of Example 4.2 for different points in the domain with different values of fractional order.For the special case of fractional order α = 1 , we compared the numerical solution obtained by YTADM with the solution in the [63] and the exact solution.The obtained result indicates that the solution obtained by YTADM has a better agreement with the exact solution than the approximate solution obtained by AVIM.One of the drawbacks of an AVIM in [63] is that the author only considered the first four iterations.It is well known that to get accurate solution considering more number of iterations is advisable.For this example, we considered the first six terms of the solution components and we obtained a more accurate solution than the one found in [63] as shown in Table 2.
Example 4.3 [50] Consider the following nonlinear time-fractional Fornberg-Whitham equation: The exact solution of this problem is u(x, t) = e x 2 − 2t 3 for the special case α = 1.
Step 1: Applying YT on both sides of (41) and using the differentiation property, we have.( 42) In our case 0 < α ≤ 1, (42) is reduced to Step 2: Applying the inverse YT on both sides of (43), we obtain.
Step 3: ADM defines the solution by an infinite series of linear terms as ( 43) and the nonlinear term as where A n , B n and C n are the Adomian polynomials to be determined.
Step 4: Using the linearity property of the YT and substituting ( 45)-( 48) into (44), we obtain (46)  Step 5: Equating the terms on both sides of (49), we have the following relation By using (9), the few components of the Adomian polynomials for the nonlinear terms uu xxx , uu x and u x u xx have been derived in the form and so on.
and so on.
Step 6: The numerical solution of YTADM is given as Here we also investigated the behavior of the approximated solution by up to three iterations of YTADM given in (50) by varying the values of the time-fractional order ( α = 0.75, 0.85, 0.95, 1 ), and the results are presented in Fig. 5, and the comparison in Fig. 6  In Table 3 we compared the obtained results of YTADM with the results of the residual power series method (RPSM) found in [50].The results show that, with error extremely close to zero, our technique provides a superior answer than the numerical method in [50].

with initial conditions
The exact solution of this problem is u(x, t) = x 2 e t for the special case α = 1.
Step 1: Applying YT on both sides of (51), we obtain.
Step 2: Applying the inverse YT on both sides of (53), to get.
Step 3: The ADM defines the solution by an infinite series of linear terms as and the nonlinear terms as where A n is the Adomian polynomial to be determined.
Step 4: Substituting ( 55) and ( 56) into (54), we obtain. ( T Step 5: Equating the terms on both sides of (57), we have the following relation.By using (9), the few components of the Adomian polynomials for the nonlinear term u 2 have been derived as Then, the iteration solution of YTADM is obtained as follows In a similar manner, we obtain.
Ŵ(3α+1) , And so on.Step 6: The numerical solution of YTADM is given as Here, we investigate the behavior graphical representations of the approximated solution of YTADM in Fig. 7,  8.The approximate solution presented in Fig. 7a and b for various values of fractional order α(α = 0.75, 0.85, 0.95, 1) is decreases as the fractional order α increases.From the comparison in Fig. 8, we conclude that the obtained results of YTADM have a nice relationship with the exact u 0 (x, t) = x 2 , A 3 (u) = 2u 0 u 3 + 2u 1 u 2 , A 4 (u) = 2u 0 u 4 + 2u 1 u 3 + u 2 2 , and so on. .
Table 4 shows the exact solution, the approximate solution of YTADM for fractional order α = 0.9 and 1 , and a comparison of the exact solution and the approximate solution obtained by HPTM, ADM, and VIM found in the reference [66] for α = 1 .The approximate solution found by YTADM is generally in considerably better agreement with the exact result than the numerical approach, in [66] as can be seen from the Table 4. Additionally, it is evident from the table that we get more accurate result as α → 1 , which indicates that the other approximate solutions are likely to be valid.

Conclusion
In this manuscript, we investigate the numerical solution of NTFPDEs using YTADM.The procedure is understandable to the readers because it consists of the direct implementation of the YT on the portion containing the fractional derivative of the given problem to change it into the algebraic form.Finally, the ADM is applied to decompose the nonlinear portion and provides a series solution to the given problem.We offered the recommended method's stability and convergence criteria along with its proof.The theoretical explanation of the suggested strategy was supported by the presentation of four illustrative instances.The findings presented in terms of figures confirm that the results obtained by the present method are in good agreement with the exact solutions for especial case α = 1 and the proposed approach gives better solution with error much near to zero.It is noted that the behavior of the approximate solution values at various α 's is identical to that of the values obtained with the exact solution, where α = 1 .This that the approximate solution is efficient.Tables 1, 2, 3 and 4 illustrate the findings, which indicate that the suggested YTADM approach outperforms the numerical methods that have been examined and published in the literature.Consequently, we draw the conclusion that YTADM is highly effective and potent in locating numerical solutions for a broad range of NTFPDEs.Only time fractional non-linear PDEs in one-dimensional spaces were used in the current work.This may be extended to the spacetime fractional PDEs with multi-dimensional spaces as well as the space fractional PDEs.Furthermore, for the purpose of addressing non-linear real-world problems, the fractional order derivative is taken in the Caputo meanings.This may also be extended to other fractional derivatives such as Caputo-Fabrizio, conformable, and Atangana-Baleanu (ABC) fractional derivative orders.

Fig. 4
Fig. 4 Solution plots of YTADM for Example 4.2 a surface of the exact solution, b surface of the approximate solution, and c surface of the absolute error, d comparison at x = 1 , e absolute error at x = 1

Fig. 8
Fig. 8 Solution plots of YTADM for Example 4.4 a) surface of the exact solution, b) surface of the approximate solution, and c) surface of the absolute error, d) exact vs approximate at x = 1 , and e) the corresponding absolute error at x = 1

at α = 1 Figure 6
shows the comparison of the numerical solution obtained by YTADM and the exact solution in the 3D and 2D plots of Example 4.3.The 3D and 2D plots have confirmed the closed contact between the YTADM and the exact solutions of Example 4.3.

Table 1
Approximate (Approx.) and comparison solutions of YTADM with the exact solution for Example 4.1 at = 0.01

Table 3
Approximate and comparison solutions of YTADM, and the exact solution for Example 4.3 when t = 0.01

Table 4
Approximate and comparison solutions of YTADM with the exact solution for Example 4.4